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Collisional damping of the excitations in a Bose-condensed 
gas is investigated over the wide range of energies and temper- 
atures. Numerical results for the damping rate are presented 
and a number of asymptotic and interpolating expressions for 
it are derived. 



The experimental realisation of Bose-Einstein conden- 
sation in dilute atomic gases [Q has led to a new range 
of experiments investigating the properties of the con- 
densate. Recently, low energy excitations of a Bose-gas 
in a magnetic trap have been studied by modulating the 
trap potential The frequencies of the lowest modes 
agreed well with theoretical predictions based on mean 
field theory [Q. Higher frequency excitations could be 
studied, for example, by light scattering Q]. Very recent 
experiments have extended the study of low-lying col- 
lective excitations to include higher temperatures where 
a finite non-condensate component may interact with the 
condensate. 

A full study of the excitation spectrum, relevant to the 
recent experimental advances, should include lifetimes of 
excitations at different temperatures and over the full 
energy range. In particular the behaviour of the damping 
rate as a function of temperature needs to be calculated. 

There have been a number of calculations of the 
condensate excitations for homogeneous systems, dat- 
ing back to the sixties and seventies [|6[-p^. These 
calculations were mainly concerned with understanding 
the physics of the phase transition. Explicit expres- 
sions where found for low-momentum asymptotics of the 
dumping rate. Its high-momentum asymptotics at T = 
were calculated by Beliaev 

In this paper, we present numerical results for the re- 
laxation of the condensate excitations which encompasses 
the whole range of energies and temperatures, and derive 
a number of asymptotic and interpolating expressions for 
the damping rate. Our results are restricted to the ho- 
mogeneous case. 

We consider Bogoliubov's quasiparticles, described by 
the field operator bp in the momentum representation. 
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1 — ttp. Here, Op is a particle field 
operator, and ap is the parameter of Bogoliubov's trans- 
formation, ap = 1 -I- p^/(2mno[/o) ~ ^p/inaUa), where 
ftp = •\/p2/2rn (p2/2r7i + 2Uono) is the quasiparticle en- 



ergy. Uq is the parameter of the collisional interaction 
of particles, written in the S'-wave scattering approxima- 
tion, Hcoii = (f/o/2) Epi-Hp2=p3+P4 4i«P2«P3ap4, ^hcre 
Uq ~ A'Ka/m, a is the scattering length and m is the mass 
of the particle. We use units where ?i = fcs = 1- 

After the Bogoliubov transformation, and neglecting 
interactions of the quasiparticles not involving the con- 
densate, we find the Hamiltonian in the form H = 
Ho + Hint, where Hq = Z^p^p^p^p and 

Hint = Uoy/n^ ^ apiOpaflpa + H. C. 

Pl+P2=P3 

= Uo^/r^; ^ ^Kn{ujp^,ujp^,ujp.^)bp^bp^bl^ 

Pl+P2=P3 

+Ka{ujp,,ujp^,u;p^)bp,bp^b-p^ + H. c. (1) 

(in all sums over quasiparticle states zero momenta are 
excluded). For brevity, we omit explicit formulae for the 
normal and anomalous interaction formfactors k„ and 
Ka . They are expressed in terms of the three scaled 
quasiparticle energies, uip^ ~ ilp^/ilo,k = 1,2,3, using 

ap = ^1 + ujp — iVp, where ilo = Uono — A-Kano/m is the 

energy scale introduced by the Bogoliubov's transforma- 
tion. 

To justify the S'-wave scattering approximation, one 
needs a^n^ < 1. Since t/onj/^ = {'ii: /mf/'^{^Qa^nQfl'^ , 
this also justifies the limit of non-interacting quasiparti- 
cles, when iJint ^ while is fixed, and allows one 
to consider Hint as a small perturbation. Note that, in 
practice, the density of the condensed phase and the tem- 
perature of the sample are functions of experimental con- 
ditions , rather than the former being a function of the 
latter. We therefore regard no as an independent param- 
eter. 

The system can be characterised by the normal average 
and the linear response function , 

Np{t~t') = lbl(t')bp{r)), (2a) 
Kp{t-t') = -ie{t)(^[bp{t):bl{t')\). (2b) 

We neglect anomalous averages (due to the anomalous 
bbb interaction) because they can be important only at 
extremely low energies. 

For non-interacting quasiparticles. 



4")(^) = -^0(^)* 



(3) 
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where rip = 1/ (e^p^^ — l ) . The chemical potential of 
the quasiparticles is zero llq] since their number is not 
an integral of the motion. K^-' (t) is the retarded Green's 
function of the 'free' Schrodinger equation, 

{zd/dt-n^)K^°\t)^Sit). (4) 

In order to derive equations for Kp and iVp, consider 
the Dyson equation for the two-point quantum averages 
in Perel-Keldish techniques , 
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Here, a, a', a", a" 
e± = Ti, 



ec^ne^n, G^""" s^""'" G^'""' . (5) 



are the C-contour indices [ll6| 



{t-t') = (TMtSbl{C)). 



(6) 



where Tc is the C-contour ordering of the field opera- 
tors and G^"^ = Gp|//,„,=n- is the exact 
self-energy. For brevity, we omit time integrations, re- 
garding the Green's functions and self-energies as kernels 
of integral operators in respect of their time arguments. 



Since 



i[Kp{t - t') - K:^{t' - t)], the 
can be expressed as linear combi- 



Green's functions (^) 
nations of Kp and A^p. The Dyson equations then 
become, 



kp = k'^^ + k^°^kpkp, 



(7) 

Np = 7V^) + k^^^kpNp + k('^apkl + ^^)«Ut (8) 



where 



Kp{t) 



(9) 



and CTp (t) is also a certain linear combination of the self- 
energy components. 

A deeper insight ||l3| shows that, ultimately, separation 
of the equation for A'p is due to microscopic causality. It 
is also very important from a more practical viewpoint. 
The initial condition Kp{0) — —i is independent of the 
interaction. Solving hence implies evolving the sys- 
tem during a finite time, and a certain simple approxima- 
tion (e.g., one- loop |l^]) to the susceptibility Kp may well 
suffice. Conversely, the knowledge of A'p implies that of 
the steady-state solution, so that simple approximations 
to dp are hopeless (cf the fact that noise sources in kinetic 
equations cannot be found perturbatively [pT|). The way 
around this problem is to find Kp and then recover Np 
using Kubo's fluctuation-dissipation theorem, thus mak- 
ing redundant. In this paper, we confine our attention 
to Kp. 

By making use of (^), Eq. |^ may be re- written in the 
Markov approximation (discussed below) as. 



( 



id/dt ~ rip + ~ Ap) Kp^S{t), (10) 



where 



(11) 



Ap-i^ = lim / dte(^">'-'^)*Kp(i). 

Taking the self-energy in the one-loop approximation px| ] 
and expressing Gp°''"" by A'p^'' and A^p"'' , after some al- 
gebra we find ||ll 7p = 7p + 7p 7p + 7p ' + 7p"- Here, 



(12a) 
(12b) 
(12c) 
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U}{LOp + Uj)k'^{uJ, OJp, UJp + U}) 



(13) 



P = \p\, Po = \/2mUonQ = y^SwanQ, 70 = Uopl/2'K = 
rJoVl287ra3no and 9 = T/VLq = 2Tm/pl. 

Two types of collision processes contribute to 7p. The 
'probe' quasiparticle can collide with a condensate par- 
ticle, producing two quasiparticles. If the thermal pop- 
ulation of the final states is neglected, this is in essence 
a classical collision, responsible for 7p. Bosonic stimula- 
tion of this process by the thermal population of the final 
states results in 7p'. The 'probe' quasiparticle can also 
collide with another quasiparticle, the final state being 
a condensate particle and a quasi-particle. This process 
is due to the bosonic attraction of the condensate, and 
results in 7p". Note that both processes contributing to 
7p are purely quantum. 

The result of direct numerical evaluation of expres- 
sions (|2|) is shown in Fig. 0. We see that both 7° and 
7p vanish if p ^ 0, hence so does 7p. This is important 
for consistency, because (i) the condensate should not be 
damped and (ii) the low-energy excitations are physically 
indistinguishable from it. Note that for low energies the 
thermal contribution always prevails. For higher ener- 
gies, 7p grows monotonically while 7p has a maximum 
at a certain momentum. It is interesting that, as the tem- 
perature grows, this maximum stabilises at fip ~ JIq, not 
at f2p T, as might be expected. The maximal value of 
7p is close to 706* = T^/1287^a^no < T. For energies high 
enough, the thermal contribution becomes negligible. 

It is easy to check that k^(ix'i, 012, cji + UJ2) ~ 

9LJiLJ2(t^l +ti^2)/32 if UJi + U2 ^ 1, ~ y2 if O^i > 1 

and w cjiy2(3 + t/|) /8(l + y|) if uji <C ^2,1, where 
j/2 = {\/^2 + 1 ^ l)/w2. We then have. 
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p<Po, (14) 



2 po 



^ ^ (1 - 5p) = 87ra2„o^ (1 - <5p) , p » po, (15) 



where (5r, 



— ^^^^^Jp^^- Except 6p, these are the well 
known results of Beliaev Q . Since the expression for 5p 
is valid only for large momenta, we are free to 'tune' it 
at low momenta so as to improve its agreement with the 

exact numerical result. On 'tuning', So ~ ^ l"!^:; {^".^ n ■ 

^ ' In 2 

In Fig. ||, we compare the results of the direct numeri- 
cal calculation with the approximate expressions. We see 
that relation (|l5|) gives a good approximation to the nu- 
merical result for p ^ po- (It even correctly reproduces 
the p^ law for low momenta.) 

Consider now the thermal contribution to the width. 
For ujp <C 1, 0, we find 7T sa 'yof{0)p/po, where 



dhJ- 



0/8 



i(cc;2 + l)(l + j/2)2^e 



-1 



(16) 



and y — (vaJ^TjTJ _ l)/a;. For low temperatures 6^1 
(T < STrano), f{e) = (37r4/20)6l4, and 



rp 37r^0^7o p inpT'^m^ 
'^P ~ 20 p'q ^ 320a2n3 



(17) 







7p dominates if Vlp 



< T. With the sound veloc- 
ity u = 2^7rano/m, ( p^ ) coincides with 7p/2 = 
'A'K^pT^ / AQriQunu^ , found by Popov [Q. For high tempera- 
tures 61 > 1, f{e) « 0.600 (cf 0) and 7^ - Qm-foOp/pa. 
Note that if ^ 1, the thermal contribution prevails for 
low momenta p <C po (cf Fig. |l|). 

A simple semi-quantitative expression for 7p, valid 
if 0,ajp » 1, can be obtained by dropping th e fac - 



t ors ff {ui^ujp) and g [uop — w, Wp) in the integrands ( |l2b[) , 
(12c), which truncate the integrals in the low-energy re- 
gion w ^ 1, setting instead the lower limit to some e > 0. 
The integration is easily performed; we then find that we 
must chose e = 2 to match the asymptotical behaviour 
of 7J at p — > 00. Then, 
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Thus the thermal contribution dominates if fJp < 
Tin (rm/87rano) T |ln na"^], where n is the total den- 
sity of particles. The last inequality is due to the fact 
that the temperature should be below the condensation 
point, T < To ~ n2/3/m. 

In Fig. ^ we compare the exact numerical results for 
the thermal component of the width with the approx- 
imate expression for low and high momenta. We have 
used a 'tuned' expression to give better accuracy at av- 
erage momenta, p ^ po, namely, 
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This expression coincides, to a very good accuracy, with 
the numerics for Wp > 1 and 9 > 10. Together, expres- 
sions (^) and ( [l9| ) provide a good approximation to 
if the temperature is not too low, 9 > 10. 

Consider now the validity of the Markov approxima- 
tion. Assume that the temperature is not too low, 
T > r^o- The bosonic distribution Hp is divergent at low 
energies. This is overcome by the 'truncating' factors 
g {uj,LUp) and g {ujp — LU^cOp) in the integrals, resulting in 
the major contribution to coming from the energies 
~ Qq. Ipso facto, it can only be applicable to time scales 
longer than l/flo- With the maximal value of 7p ~ 706*, 
our results apply if ^ Qq/^q. This can be written as 



(20) 



To understand what this restriction means in terms of 
the critical temperature Tc, note that rnT ~ n'2/^, where 
n' is the density of the non-condensate phase. Estimate 
( pO| ) then is equivalent to n' <C no/{nQa^)^^^. This does 
not contradict n' ~ rig. Thus, although certainly T < Tj,, 
a stronger condition T <^Tc does not seem necessary. 
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FIG. 1. Results of the direct numerical calculation of the 
width components 7p (dashed line) and 7p (sohd lines). 

FIG. 2. Comparison of the results of the direct numerical 
calculation of 7p (solid line) with the approximate expres- 
sions: dashed lines - Beliaev's asymptotics for small and large 
momenta, dash-dotted line - 'tuned' expression (|l5|). 

FIG. 3. Comparison of the results of the direct numerical 
calculation of 7p (solid lines) with approximate expressions 
(^^ (dashed lines) and ( |l9| ) (dash-dotted lines). 
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